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Abstract. A possible slowing down of the cosmic expansion is investigated through 

^^ . a cosmographic approach. By expanding the luminosity distance to fourth order and 

r) I fitting the SN la data from the most recent compilations (Union, Constitution and 

^ ' Union 2), the marginal likelihood distributions for the deceleration parameter today 

Oj. suggest a recent reduction of the cosmic acceleration and indicate that there is a 



considerable probability for qo > 0. Also in contrast to the prediction of the ACDM 



\—( • model, the cosmographic q{z) reconstruction permits a cosmic expansion history where 

C/3 , the cosmic acceleration could already have peaked and be presently slowing down, 

^ ' which would imply that the recent accelerated expansion of the Universe is a transient 



phenomenon. It is also shown that to describe a transient acceleration the luminosity 
distance needs to be expanded at least to fourth order. The present cosmographic 



> 

vQ ' results depend neither on the validity of general relativity nor on the matter-energy 

OO . contents of the Universe. 

a^ 

l/^ ■ Keywords: supernova type la - standard candles, dark energy experiments, dark energy 

S ■ theory 



1. Introduction 

Independent measurements by various groups have suggested that the recent expansion 
of the Universe is speeding up and not slowing down, as believed for many decades 
[Ij. In other words, by virtue of some unknown mechanism, the cosmic expansion 
underwent a "dynamic phase transition" whose main effect was to change the sign 
of the universal deceleration parameter q{z) at low redshifts. Together with other 
complementary observations, such results are strongly indicating that we live in a flat, 
accelerating Universe composed of ~ 1/3 of matter (baryonic -|- dark) and ~ 2/3 of an 
exotic component with large negative pressure, usually named dark energy (see [2] for 
some reviews). 

However, what seems more certain is that the Universe was accelerating in a recent 
past, but that its current state of acceleration is less well determined [3]. Some authors, 
using a dynamical ansatz for the dark energy equation of state, have suggested that the 
cosmic acceleration has already peaked and that we are currently witnessing its slowing 
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down [1]. This behavior was also suggested by other works combining different methods 
and data f5]. 

The possibihty of a transient cosmic acceleration appears in several theoretical 
scenarios [6] and is also theoretically interesting because eternally accelerating universes 
(like ACDM) are endowed with a cosmological event horizon which prevents the 
construction of a conventional S-matrix describing particle interactions. Indeed, such 
a difficulty has been pointed out as a severe theoretical problem for any eternally 
accelerating Universe [3 E]. Naturally, whether the acceleration of the Universe is 
slowing down the overall dynamic behavior in the future may be profoundly different 
from the one predicted by the ACDM evolution. 

Here we investigate the question related to a possible slowing down of the recent 
accelerating stage of the Universe by using a distinct method. In order to access the 
present value of the deceleration parameter we consider the so-called cosmographic 
approach [9l [TOl [HI [121 [13 [IH [13 [IS] instead of assuming the validity of a specific 
gravitational theory together with a given dynamical dark energy model. As happens 
with kinematic models [3l[Tl[I71[l8l[l9l[20l[2ll|22l[23l[2l], the basic advantage of 
the cosmographic approach over dynamic models is that the former depends neither on 
the validity of general relativity nor on the matter-energy contents of the Universe. In 
addition, to verify the dependence of the method on a particular set of supernova data, 
we repeat our analysis by considering three different samples, namely, the Union |25] . 
Constitution [26], and Union 2 [27] compilations. 

2. Cosmographic approach 

Following a traditional approach, the late time cosmic expansion can be expanded as 
[Tni[28] 



a{t) = l + Ho{t-to)-^qoH^oit-tof + l^JoH^it-tof + ^soH^{t-to)^ + 0[{t-to)%{l) 

where Hq, qo, jo and sq are the present day values of the Hubble, deceleration, jerk and 
snap parameters, respectively. Similarly, the luminosity distance can also be expanded, 
yielding an extended version of the Hubble law 



dL{z) = — 
-no 



z + -(1 - qo)z^ - -(1 - go - 3go' + Jo)^:' 



+—{2 - 2qo - 15ql - 15q^o + lOgoJo + 5jo + So)z' 



+ 0(/), (2) 



where it was assumed a flat universe, as is done throughout this work. As it appears, 
the truncation of the above expression is endowed with some convergence issues [12] . 
which is a source of concern since the SN la samples do have events at ^ > 1. One 
way to circumvent this problem of convergence at high redshift is to limit the supernova 
sample to events of lower redshift. Another approach is to parametrize cosmological 
distances by the so-called y-redshift [T2], defined hj y = (Aq — Ae)/Ao, a quantity related 
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to the usual redshift, z = (Aq — \e)/\e, by 

Using the y-redshift we can expand the luminosity distance as 



ddy) = TT 

-no 



1/ + 2 (3 - 9o)y^ + g(ll - 5go + 3go' - Jo)r 



+ ^(50 - 26go + 21^0^ - I5g3 + iQg^j^ _ 7^^ + so)y' 



+ Oiy'U) 



This parametrization has the nice property that its convergence radius comprises the 
full range of past cosmological events. 

In what follows, in order to fit the SN la data compiled by the Union (307 events), 
Constitution (397 events) and Union 2 (557 events) samples, we consider (i) the fourth 
order expansion of the luminosity distance in redshift, which we will label d\^ (z), that 
is, the truncation of ([2]), and (ii) the fourth order expansion of the luminosity distance 
in y-redshift, d\^ (y), the truncation of (jll). For the redshift expansion fit we use sub- 
samples with z < Zcuti labeling the expansions d)^ {z < Zcut), which for Zcut = 0.5 yields 
157 events (51%) for the Union, 247 events (62%) for the Constitution, and 402 events 
(72%) for the Union 2 compilation. 

We construct the likelihood as £ = exp[— x^/2], where 

2 _ V^ [fiBii;MB,a,l3) - fi{i; Ho, qo, jo, So)f 



X 



^^77^ ' (^) 



^;i« 



where i runs through the SN la events, /x = 51ogiQ[(ii/(lMpc)] + 25 is the predicted 
distance modulus in a given model, and o"^ is the total uncertainty on the distance 
moduli. In the likelihood analysis we analytically marginalize over [211 129] the 
nuisance parameter containing the Hubble parameter (Ho) and the supernova absolute 
magnitude {Mb), and numerically maximize for two remaining nuisance parameters 
related to the supernova light curve fitting {a and /3). For both the Union and 
Constitution sets the SALT method values for m^°-^ , s (stretch) and c (color) were 
used: hb = mJ^'^^ — Mb + a(s — 1) — /3c. Union 2 uses SALT2 light curve fitting: 
Hb = m^°'^ — Mb + cxxi — /3c, where xi is a free parameter that measures the deviation 
from the average decline rate. 

3. Current deceleration probability 

Since we are mainly interested in the value of the deceleration parameter today, we 
show in figure [H the marginal likelihoods for qq (marginalization over jo and sq). Both 
cosmographic parametrizations, d\^ [z < 0.5) |j and d\^ {y), are statistically compatible 
and have similar dispersion, even though the 2 < 0.5 implies a considerably smaller 
number of events in the d\^ {z < 0.5) analysis in relation to the full SN la sets used in the 

I The choice of Zcut = 0.5 here is meant to be conservative in terms of possible convergence problems 
in the d]^ {z) expansion. 
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Figure 1. Marginal likelitiood for the deceleration parameter today. Top panel uses 
the Union SN la compilation, central panel uses the Constitution, and the lower panel 
uses the Union 2. In the upper left corner of each panel is indicated the number of 
events in each data set. Curves are scaled to ina.x[C{qo)] = 1. 



d^ {y) analysis. This effect was already noticed before [13], and that may be indicating 
that the y-redshift expansion is not optimal to maximize the Fisher information. 

The basic distinction feature between the Union and Constitution sets comes from 
90 events at low redshift {z < 0.5). They have the effect of shifting the go likelihoods 
to higher values (figure [1]) , thereby increasing the probability of a decelerating universe 



Could the cosmic acceleration be transient? A cosmographic evaluation 
Table 1. Probability of a positive deceleration today, P{qa > 0). 



fit \ set 


Union 


Constitution 


Union 2 


d^\z)l 


0.006 


0.097 


0.008 


df{z < 0.9) 


- 


- 


0.082 


df{z<Q.l) 


- 


- 


0.40 


d^\z< 0.5) 


0.08 


0.27 


0.55 


di\z< 0.3) 


- 


- 


0.57 


df{y) 


0.15 


0.65 


0.33 



(table [I]). The Union 2 has an even higher percentage of low redshift events and the 
deceleration likelihood derived from d\^ [z < 0.5) is further shifted to higher values of 
go- This suggests that the low-z SN la events are pointing to a higher go than the full 
sample. Such a trend could help to explain why Shafieloo et al. |1] found a slowing 
down of the cosmic acceleration when analyzing the Constitution data set, but probably 
would find (i) a weaker result had they used the first Union compilation and (ii) maybe 
an even stronger result had they used the Union 2. 

To further explore the suggestion that low-z SN la events may be pointing to higher 
go values, we examine the Union 2 data set up to various redshift cuts, Zcuti with the 
dj^ {z < Zcut) parametrization. This point of view is reinforced by the results shown 
in figure [2] and table [H The more recent the SN la events are, the higher is the peak 
and average deceleration values today and the larger is the probability for positive go. 
Note that the higher deceleration probability with lower redshift cutoff is due not only 
to a higher peak and average deceleration parameter, but also to a higher variance in 
the distribution, consequence of looser constraining of the parameter space. The use 
of various redshift cuts in the d\^ (y) parametrization also corroborates those findings 
(numerical results not shown). 

In table [H by using d\^ {z < Zcut) and d^^'{y), we quantify the probability for the 
universe to be slowing down today for the three SN la samples. It has been estimated 
by computing the integrated marginal likelihood for go > 

When just z < 0.5 events are fitted, then the Constitution and Union 2 sets indicate a 
considerable probability of a slowing down universe. This probability is also large (and 
even dominant in the Constitution case) when the samples are fitted by d\^ (y). Using 
a fifth order expansion, John [9l[T6] also found a significant probability of a positive go. 
It is not totally surprising that both parametrizations, d}^ {z < Zcut) and d\^ {y), 
lead to distributions for go that are not identical. Previous works (e.g. [Hj) showed 

§ The use of z > 1 SNe in the d\^ {z) parametrization incurs more severe convergence problems, 
therefore these results must be taken with greater caution. 
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Figure 2. Probability distribution function for tlie deceleration parameter today for 
the d]^ {z < Zcut) parametrization. The internal legend shows the cut redshift and the 
corresponding number of SN la from the Union 2 compilation. 



that best-fit values are model (parametrization) dependent, even for models with same 
goodness of fit. 

Because the d^ {z < 0.5) and d\^ [y) fits yield broad marginal likelihoods that are 
nearly symmetrical and peaked around go = O5 most notably with the Constitution and 
Union 2 compilations, P(go > 0) is very sensitive to the point of maximum likelihood, 
which is not a robust feature of those fits in relation to the various data sets. This 
explains the variability of the results shown in table [1] in relation to the different SN la 
sets. 

As a byproduct, we also obtain the marginal likelihood distributions for the jerk 
and snap parameters today (figure [3]) from the Union 2 data. The dj^ {z < 0.5) 
parametrization yields somewhat tighter constrains than dj^ (y), even though it uses 
less events. Our results for the y-redshift parametrization are also in general agreement 
with what [15] found using the first Union data set. For reference, the flat ACDM model 
has go = |^m — 1, jo = 1 and Sq = 1 — ^fim- Therefore this model has a single free 
parameter, which the best fit to the Union 2 data gives go = —0.565 ± 0.032. 



4. Cosmographic deceleration history reconstruction 



Independently of the cosmic deceleration today, it is of interest to know if there is another 
change in the sign of the cosmic acceleration after a prior transition from a decelerated 
to an accelerated phase at moderate redshifts {zt ~ 0.5 — 1). From the power series 
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Figure 3. Marginal likeliliood for tire jerk and snap parameters today from the Union 
2 data. 



expansion of the scale factor ([T]) one can also express the deceleration parameter as a 
power series in time. This time variable can be written as a power series in redshift or 
y-redshift, yielding respectively 



q{z) = go + (-go 
and 

Qiy) = qo + 



2ql+jo)z + -{2qo 



1 



^^0 + 8^0 - 7goJo - 4jo - so)z^ + 0{z^ 



go - 2go + jo)y + -(4go + Sq^ - 7gojo 



2jo - so)y^ + 0{y' 



(7) 



(8) 



Equivalently, if the truncation of ([T]) is made at third order, then the resulting expression 
for q{z) and q{y) is linear. Several works [HI [ill EOl EH |24] were based on the 
q{z) = qo + qiz parametrization for the deceleration, which does not allow a transient 
cosmic acceleration. The derivation of a power expansion for the deceleration from the 
scale factor expansion establishes a natural link with the kinematic studies based on 
a q{z) expansion in power series. J. V. Cunha (private communication) employed the 
q{z) = go + qiz + ^g2-2^^ parametrization and found a suggestion of transient acceleration 
(see [22] for another formulation of a quadratic kinematic model). The quadratic form 
of ([7]) and ([8]) allows for a decelerated past, a transition to an accelerated phase, a point 
of maximum acceleration, then a slowing down of the acceleration and a transition 
to a recent or future decelerating phase. Similar behavior of transient acceleration is 
predicted or allowed by several dynamic models [HE]- In contrast, the ACDM model 
predicts a monotonic deceleration history connecting its asymptotic limits in the past 
and future, q{z — )■ oo) = 0.5 and q{z -^ —1) = —1. 

Even though ([7]) and ([8]) are expansions of same order in their respective variables, 
the second is a more precise expression since the order of its error 0{y^) = 0{z^/[l + z]^) 
is always lower than the order of the error for the first, 0{z^), in the past. 

We use the cosmographic parameters (go, jo and Sq) obtained in each 
parametrization with the Union 2 sample to reconstruct q{z) using (jH]) and ([3]), which 
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Figure 4. Deceleration liistory reconstruction. The gray area and the darker area 
inside it delimit the la confidence region for the q{z) reconstruction obtained from 



various Union 2 SN la analyzes. Top panel: d\^'{z < 0.3) (gray), d\^'{z < 0.7) (darker 
internal area) and flat ACDM model (green thinner area). Bottom panel: d]^ (y) with 
no extra conditions (gray) and with the condition that the deceleration parameter must 
be positive at z = 2 inside the la confidence region (darker internal area). 



is shown in figure HI The confidence regions shown in the figure are drawn from all 
q{z) curves allowed by the parameter values such that C{qQ, Jq, Sq) / Cmax is less than a 
threshold [30] . The top panel shows the deceleration reconstruction for the d^ {z < Zcut) 
parametrization for Zcut = 0.3 and Zcut = 0.7. The first constrains very loosely q{z), not 
even showing a phase of accelerated cosmic expansion with confidence above la. The 
higher cut redshift allows a tighter constraining, being maximum for 2; ~ 0.1. 

The q{z) reconstruction obtained from the d\^ {y) parametrization is shown in the 
bottom panel of figure HI It defines a clear phase of accelerated expansion around 
z ~ 0.2, but suggests it is slowing down around z ~ 0.1, where the q{z) reconstruction 
is maximally constrained. 

For comparison, figure H] also plots the results for a fiat ACDM model, showing 
that under this parametrization all redshifts have equally well determined deceleration. 
As we know that the SN la data is limited to a redshift range, it is clear that the 
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preponderant factor in this ACDM result is the model itself, not the data. We also 
note that, even though the comparison of the cosmographic parametrizations with the 
ACDM model is not the objective of this work, all have very similar goodness of fit as 
measured by the reduced x^. 

All cosmographic parametrizations contain transient acceleration solutions in their 
la confidence regions with a considerable positive deceleration probability today and 
in the near future, if we assume that the cosmographic parametrization has some 
predictive power. However, due to the large uncertainties in the q{z) reconstructions 
at large redshifts or in a projection into the future {z < 0) from the cosmographic 
parametrizations of the Union 2 data alone, it is not possible to affirm if there is 
a transition from a phase of accelerated expansion to a present or future phase of 
deceleration. 

An improvement in the containment of the allowed q{z) could be achieved by 
imposing external conditions to the cosmographic analysis, for example that in the 
remote past the matter dominated universe should be decelerated. As an exploratory 
investigation, we implement this external restriction to the analysis in the bottom panel 
of figure m by imposing that q{z = 2) > 0[||1 at the la confidence level. As a result, the 
allowed q{z) is further constrained, only leaving transient acceleration solutions. 

The strict transient behavior obtained with the external q{z = 2) > condition may 
be a consequence of the fourth order cosmographic parametrization for the luminosity 
distance, as this implies a quadratic q{z) which may be forced to have two zeros (two 
transition redshifts). Every parametrization implies a bias, and to test if this is the case 
for the transient solutions in the approach adopted, one would need to go to higher order 
in the cosmographic expansion. However, models with a higher number of parameters 
would need more data to meaningfully constrain the parameter space. 

5. Conclusion 

To address the main question of this work, we should focus on the results for the 
Union 2 compilation, since it is currently the largest and most updated. Both kinds 
of cosmographic parametrization, d\^ {z < z^ut) and d]^ (y), allow the circumvention of 
convergence issues related to the power expansion, broadly agree and (i) suggest that 
lower redshift SN events point to lower cosmic acceleration, and (ii) indicate that there 
is a considerable probability for the Universe to be decelerating today. The slowing 
down of the cosmic acceleration and even a transient acceleration scenario is allowed in 
all cosmographic parametrizations. 

The possibihty that those results may be due to bias induced by the form of the 
parametrizations and data analysis artifacts may be examined inside the cosmographic 
approach by going to higher order in the scale factor expansion and through the use of 

II The choice of z = 2 is arbitrary, but meant to be a conservative comproniise. A larger redshift would 
imply a lower constraint on the fitting, and a with a smaller redshift, the constraint would be stronger, 
but then we would be less certain about the reality of the q{z) > condition. 
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simulated data inside some conveniently chosen fiducial cosmological models. That in 
itself would constitute an interesting follow up study that could be useful to improve 
the understanding of the limitations of the observational data sets and of the analysis 
methods. 

What may be causing a possible slowing down of cosmic acceleration and transient 
cannot be answered with basis only in a cosmographic approach. We recall that in 
the accepted scenario (inflation followed by radiation, dark matter and dark energy 
dominated stages), the Universe already had two accelerating stages and what is being 
discussed is whether the present supernova data are compatible with a new transition to 
a decelerating phase. Some examples of dynamical scenarios discussed in the literature 
include quintessence models in which the equation of state parameter becomes positive 
in the future, and some D-branes inspired scenarios, among others [6l|8]. Naturally, we 
cannot exclude the possibility that such a slowing down might be provoked by some 
exotic kinematic (or dynamic effect) in the presence of inhomogeneities [3T] . 

In conclusion, we used fourth order expansions of the luminosity distance to describe 
the recent cosmic expansion as probed by SN la observations. We showed that one 
needs to go at least to fourth order in the expansion of the luminosity distance to 
be able to describe a transient acceleration. Other works [131 US]) using smaller data 
sets than some of those used here, indicate that going to higher orders does not add 
meaningful statistical information. Under a cosmographic analysis of the most recent 
SN la compilations, we have found a suggestion of cosmic acceleration slowing down 
and the existence of a considerable probability in the relevant parameter space that the 
Universe is already in a decelerating expansion regime. This result is in great contrast 
to the standard ACDM model, which predicts that the cosmic expansion is accelerating 
and must remain so forever with an ever increasing acceleration. 

The predictability of the cosmographic parametrizations are considerably inferior 
to the ACDM model and other dynamical models, but that seems to be the price 
for avoiding the usual assumptions about the physical contents of the Universe and 
the underlying gravity theory. It should also be stressed that the price for avoiding 
convergence issues for the expansion in powers of z has severe consequences on model 
predictability, but this could be overcome by a much larger SN la sample (which is 
foreseen in future surveys), by the use of other kinds of data and priors, and maybe by 
a smarter (optimal) cosmographic modeling. 
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